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Abstract. We use an analytic approach to study the susceptibility of the Planck Low Frequency Instrument 
radiometers to various systematic effects. We examine the effects of fluctuations in amplifier gain, in amplifier 
noise temperature and in the reference load temperature. We also study the effect of imperfect gain modulation, 
non-ideal matching of radiometer parameters, imperfect isolation in the two legs of the radiometer and back-end 
1// noise. We find that with proper gain modulation 1/f gain fluctuations are suppressed, leaving fluctuations 
in amplifier noise temperature as the main source of 1// noise. We estimate that with a gain modulation factor 
within ±1% of its ideal value the overall 1/f knee frequency will be relatively small (< 0.1Hz). 
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1. Introduction 

■ 

^ [ PlanckQ is a European Space Agency (ESA) satellite 
' mission to map spatial anisotropy and polarization in the 
?3 , Cosmic Microwave Background (CMB) over a wide range 
of frequencies with an unprecedented combination of sen- 
sitivity, angular resolution, and sky coverage (Bersanelli et 
al. fL996p. Following the breakthrough of the COBE dis- 

Bcnnctt et 



1992 



covery of CMB anisotropy (Smoot et al. 
al. |1996p , and the MAPQ satellite launched in June 2001, 
Planck will be the third generation space mission dedi- 
cated to CMB observations. 

The data gathered by these missions will revolutionise 
modern cosmology by a precise determination of the 
fundamental cosmological parameters which govern the 
present expansion rate of the universe, the average density 
of the universe, the amount of dark matter, and the nature 
of the seed fluctuations from which all structures in the 

19951 , 



1997, Scott et al. 



universe arose (see, e.g., Hu et al. 

White et al. 1994 for recent reviews on CMB anisotropy) 
and will provide at the same time full sky surveys at es- 
sentially unexplored frequencies, with fundamental impli- 
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Planck homepa ge: http://astro.cstec.esa.nl/Planck/ 
2 MAP homepage: http://map.gsfc.nasa.gov 



cations for a large area of problems in astrophysics (De 
Zotti et al. 



199S) 



Planck consists of a High Frequency Instrument 



(HFI) (Puget et al. 1998) and a Low Frequency Instrument 
(LFI) (Mandolesi et al. 199§| ) observing the sky through a 



common telescope. While the MAP radiometers measure 
temperature differences between two widely separated re- 
gions of the sky through a pair of symmetric back-to-back 
telescopes, the Planck LFI radiometers are designed to 
measure differences between the sky signal and a stable 
internal cryogenic reference load. The LFI scheme takes 
advantage of the presence in the focal plane of the HFI 
front end unit, which is cooled to approximately 4 K as 
an intermediate cryo stage for the 0.1 K bolometer detec- 
tors. 

The LFI radiometer design is a modified correlation 
receiver (Blum 1959, Colvin 1961, Bersanelli et al. 1995), 
realised with High Electron Mobility Transistor (HEMT) 
amplifiers at 30, 44, 70 and 100 GHz. The modification is 
that the temperature of the reference load can be made 
significantly different from the sky temperature. To com- 
pensate for the offset (a few K in nominal conditions), a 
gain modulation factor, r, is used to null the output sig- 
nal in order to minimise sensitivity to RF gain fluctuations 
and achieve the lowest white and 1/f noise in the output. 
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Obtaining data streams with low 1/f noise is of pri- 
mary importance in order to achieve the LFI scientific 
objectives. In fact excessive 1// noise would degrade the 
quality of the measured data (Janssen et al. 1996) by in- 
creasing the effective rms noise and the uncertainty in the 
power spectrum at low multipole values. Such effects can 
be avoided if the post detection knee frequency fk (i.e. 
the frequency at which the 1/f contribution and the ideal 
white noise contribution are equal) is significantly lower 
than the spacecraft rotation frequency (/spin ~ 0.01 7 Hz). 
For values of fk greater than / sp m it is possible to mitigate 
such effects by applying appropriate destriping and map 
m aking algorithms^ to the time ordered data (Maino et 
al. |2000D . 



If the knee frequency is sufficiently low (i.e. fk < 
0.1 Hz), with the application of such algorithms it is pos- 
sible to maintain the increase in rms noise within few % 
of the white noise, and the power increase at low multi- 
pole values (i.e. I < 200) at a very low level (~ two order 
of magnitude less than the CMB power). If, on the other 
hand, the knee frequency is high (i.e. ^> 0.1 Hz) then even 
after destriping the degradation of the final sensitivity is 
of several tenths of % and the excess power at low multi- 
pole values is significant (up to the same order of the CMB 
power for ~ 10 Hz, Bersanelli et al., 2002). Therefore, 



careful attention to instrument design, analysis, and test 
is essential in order to achieve a low 1/f noise knee fre- 
quency. 

In this paper, we analyse the most important system- 
atic effects due to non-ideal behaviour of components in 
the LFI radiometer signal chain, and estimate their im- 
pact on the post detection knee frequency. In section || we 
present a general analytical description of the Planck- 
LFI radiometers to derive formulas for the radiometer 
power output and sensitivity in the two cases of perfectly 
balanced and slightly unbalanced radiometer. Here we also 
show that under quite general assumptions, the radiome- 
ter sensitivity does not depend on the reference load tem- 
perature. In section |^ we analyse the impact of various 
systematic effects on the post-detection knee frequency 
showing that with proper gain modulation it is possible 
to keep the radiometer 1/f noise to a very low level also 
in the presence of different non-ideal behaviours (e.g., gain 
and noise temperature imbalance, imperfect gain modu- 
lation). For sake of conciseness, we transfer part of the 
formalism to the appendices. Finally, in section ^| we sum- 
marise our results and discuss briefly their implications 
for Planck observations. 



2. Analytic model of LFI pseudo-correlation 
radiometers 



sec Burigana et al. f 997 , Delabrouille 199£ , Main o et al. 

Natoli 



199S , for d etails about destriping and Dore et al. |200l| 



2.1. Radiometer architecture 

In Figure [j] we show a schematic of the baseline LFI ra- 
diometer design. In this design, each feed-horn is con- 
nected to a Radiometer Chain Assembly consisting of an 
actively-cooled 20 K front-end connected to a 300 K back- 
end via waveguides. 



Reference @ 4K 



Radiometer front-end 



Radiometer front-end 



Reference @ 4K 



Waveguides 



Radiometer back-end 



Radiometer back-end 



et al. 2001, for details about map-making algorithms. 



Fig. 1. Schematic of the Planck LFI radiometer design. 
Each feed horn is connected to two radiometers consti- 
tuted of a 20 K front-end and a warm back end at 300 K. 
In our design the reference load has a temperature of 4 K. 

In the front-end part (see top part of Fig. ||) the radia- 
tion entering the feed-horn is separated by an OrthoMode 
Transducer (OMT) into two perpendicular linearly po- 
larised components that propagate independently through 
two parallel radiometers. In each radiometer, the sky sig- 
nal and the signal from a stable reference load at ~4 K 
are coupled to cryogenic low-noise High Electron Mobility 
Transistor (HEMT) amplifiers via a 180° hybrid. One of 
the two signals then runs through a switch that applies 
a phase shift which oscillates between and it at a fre- 
quency of 4096 Hz. A second phase switch will be present 
for symmetry on the second radiometer leg; this switch 
will not introduce any phase shift in the propagating sig- 
nal. Therefore it will not be considered in our analysis. 
The signals are then recombined by a second 180° hybrid 
coupler, producing an output which is a sequence of sig- 
nals alternating at twice the phase switch frequency. 

In the back-end of each radiometer (see bottom part 
of Fig. |) the RF signals are further amplified, filtered 
by a low-pass filter and then detected. After detection the 
sky and reference load signals are integrated, digitised and 
then differenced after multiplication of the reference load 
signal by a so-called gain modulation factor, r, which has 
the function to make the sky-load difference as close as 
possible to zero. 

According to this architecture each radiometer will 
produce two independent streams of sky-load differences; 
the final measurement is provided by a further averaging of 
these differenced data samples between the two radiome- 
ter legs. 

The LFI pseudo-correlation design offers two main ad- 
vantages: the first is that the radiometer sensitivity does 
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— x(t) and represent the voltages at the sky and 
reference horns respectively; 

— gp x , gp 2 , g Bl and g B2 are the voltage gains of front-end 
and back-end amplifiers; 

— np 1 (t), np 2 (t), n Bl {t) and n B2 (t) represent the noise 
voltages of the front-end and back-end amplifiers; 

— 4*Fn t pF 2 i <pBi an d 4>b 2 represent the signal phases after 
the amplifiers; 

— Q\ and 62 represent the phase shifts in the two switch 
states (in LFI baseline 9\ = and 62 = 7r); 

— Ai and A2 represent the fraction of the signal ampli- 
tudes that is transmitted after the phase switch in the 
two switch states. For lossless switches A\ 2 = 1. 



The output signal at the two radiometer legs is given 



by: 



Fig. 2. Details of the front-end (upper figure) and back- 
end (lower figure) of LFI radiometers. 

not depend (to first order) on the level of the reference 
signal; the second is provided by the fast switching that 
reduces the impact of 1 // fluctuations of back-end ampli- 
fiers. In the next subsection we present the analytical de- 
scription of this radiometer design. In particular we derive 
formulas for the power output and sensitivity in the two 
approximations of (i) perfectly balanced and (ii) slightly 
unbalanced radiometer. Then we show how the radiometer 
sensitivity is not dependent, to first order, on the reference 
load temperature and we also provide an estimate of this 
dependence in the case of slightly unbalanced radiometer. 

In order to improve readability we have kept the math- 
ematical treatment as simple as possible, relegating the 
full definition of long formulas to appendices. 

2.2. Analytic model of LFI radiometers 

2.2.1. Input-output signal sequence 

Referring to the schematic in Fig. |^, we can define the 
following transfer functions for the different radiometer 
components. 



/hyb : {x, y} 



/amp : {. x i U} 



x + y x — y 



{g Fl (x + n Fl )e^i , 9 p 2 (y + np 2 )e 4 ^ } 
/sw : {x,y}^ {x,y/A~e^y} with j = 1,2 (1) 

/amp • { x i V\ * 



{g Bl (x + n Bl )e^ , g B2 (y + n B2 )e 1 ^ } 



Sj,i = e^^gs, < n Bl 



V2 



e l ^g Fl [n Fl 



x + y 
V2 , 



+ X /A~e l(ej+0F 2 ) 



g F2 n F2 



x-y 

V2 



(2) 



^/2 



fA^^gp, { n,, 2 ^ 



x-y 



where the subscript j = 1,2 indicates the two phase switch 
states. 



2.2.2. Radiometer power output and sensitivity 

Our first step is to write an expression for the radiometer 
power output at a given time t, considering that in our 
scheme we take the difference between the sky and load 
signals. If we indicate with pi(t) the power output at each 
of the two radiometer legs (where 1 — 1,2) and with p (t) 
the radiometer power output we have that: 



pi(t) = a (\s 1A \ 2 (t) - r\s 2 ,i\ 2 (t + At sw )) 

p 2 {t) = a (|s 2 , 2 | 2 (i + A^w) - r\s h2 \ 2 (t)) 
p(t) = l/2[ Pl (t)+ P2 (t)} 



(3) 



where a is the constant of proportionality of the square- 
law detectors, At sw ~ 122 fis corresponds to half of the 
phase switch period and r is the gain modulation factor 
which in LFI scheme is applied in software. The expanded 



where: 



form of pi{t) is reported in equation (A.l); a similar rela- 
tionship holds for p 2 (t). The assumption of identical val- 
ues of a in the above equation can be made without loss 
of generality; any difference in actual values between the 
two legs could be absorbed into back-end gain differences 
for the purposes of this analysis. 

In our treatment we will show how the gain modula- 
tion factor (a real number generally < 1) can be tuned 
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properly in order to have an approximately null power 
output, which leads to almost complete suppression of the 
dependency of the radiometer sensitivity on the reference 
load signal, and minimal impact of gain fluctuations in 
the front-end amplifiers. Note that in equation ([|) the sky 
and load signals appearing in the difference are sampled at 
slightly different times, which are relative to the two differ- 
ent phase switch states; this is relevant when we consider 
the effect of 1// fluctuations of the back-end amplifiers. 
The effect of these instabilities can be practically elimi- 
nated by using a phase switch frequency / sw m uch greater 
than the 1// noise knee frequency; in section 3J5 we will 
show that the baseline phase switch frequency of ~ 4 KHz 
considered for LFI leads to a very low level of back-end 
f // noise. 

Let us now calculate the time average of pi(t) inte- 
grated over the bandwidth, (3, i.e. p l — L J At pi(t)dtdv. 
Considering that signal and amplifier noises are uncorre- 
cted we have that cross-correlation terms vanish; there- 
fore we can write p ± and p 2 in the following compact form: 



p 1 = akfiGB! x 



1 - 



T X (G -rI)-rT y \G--I\ + T ni - rT ni 



(4) 



p 2 = ak(3Gs 2 x 



T X (G -rI)-rT y [G--I)+ T„ 2 - rT n 



where: 



G =4 

7=1 

4 



G Fl + A 1 G F2 + 2y/ A 1 G Fl G F2 cos(0i 



G Fl + A 2 G F2 + 2y/ A 2 G Fl G F2 cos((9 2 



G Fl + A 2 G F2 - 2y/A 2 G Fl G F2 cos( 



G Fl + A l Gp 2 - 2 y /A 1 G Fl G F2 cos(9 1 



Tm — -^{G Fl T nFi + AiG F2 T nF2 + 2T nBi ) 
Tn, = -^{G Fl T nFi + A 2 Gp 2 T nF<2 +2T„ B; ) 

In equation (§) the terms T x ,T y , T nFi , T nF2 , T nB2 and 
T nB2 represent the sky, load and amplifier noise tempera- 
tures which are defined by relationships like k(3T x = x 2 (k 
is the Boltzmann constant), G Fl ,p 2t B 1 ,B 2 — 9 F f b b 
are the amplifier power gains, (3 is the bandwidth and 
4> = tf>p 2 — (j) Fl is the phase mismatch of the front-end 
amplifiers. 

In the case in which radiometer parameters (gains, 
noise temperatures, etc.) depend on the frequency, equa- 
tion (Q) is still valid provided that we use values aver- 
aged over the bandwidth, i.e. for each parameter P on 



the left-hand side of the above definitions: P — ► P = 

We now see that in order to null the output at each 
radiometer branch output (i.e. make pi — 0), we must 
adjust r to a proper value, denoted r* . In the general case 
discussed above we have that: 



P2 =0 



T X G 



r = r, = 



Tyl 



TyG 

T X G 



TJ. + T n 



Tyl 



TyG- 



(5) 



It should be noted that both the sky temperature 
and the instrumental parameters arc expected to undergo 
changes of order of few mK. Thermal variations and ag- 
ing are likely to drive instrumental changes which are ex- 
pected to occur on time scales of days to months. The main 
source of change in T x is the CMB dipole signature, which 
introduces a spin-synchronous (1 rpm) temperature mod- 
ulation ST X < 3.5 mK. In principle it could be possible to 
implement feedback schemes in order to vary r almost in 
real time and have a constant zero output , however, this 
would complicate the scheme and possibly introduce addi- 
tional sources of systematics. We anticipate that adjusting 
r on the timescale of a few days will be sufficient, so that 
we can consider r constant in the following calculations. 
This issue is discussed in more detail in a forthcoming 
paper currently in preparation. 

The radiometer sensitivity AT (i.e. the minimum sig- 
nal that can be detected in a bandwidth [3 and integration 
time r) can be calculated following the approach outlined 
in appendix [b| where we also report the complete algebraic 
form. 

2.2.3. Approximations 

In this section we derive some approximations of the gen- 
eral equations presented in appendices ^ and [b|. 

Perfectly balanced radiometer. The first, zero-order ap- 
proximation is relative to the ideal case of a perfectly bal- 
anced radiometer. We assume that the radiometer compo- 
nents are ideally matched, but that there still is a temper- 
ature offset between the sky and reference load, and that 
there are still noise and gain fluctuations present in the 
amplifiers. This approximation can be derived by setting 
G Fl = Gp 2 = Gf, Gbx — Gb 2 = Gb, T nFi = T„ F2 = T„ F , 
T nBl = T nB2 = T„ B , 6 l = 0, 9 2 = 7T, <t> = 0, A 1 = A 2 = 1. 
In this case the radiometer power output and sensitivity 
can be written in the following simple form: 



P^Po = akf3G F G B [T x + T ntot - r(T y + T ntot )] 



4 In this case the astrophysical information would be recov- 
ered by recording the gain change needed to maintain the ra- 
diometer balance. 
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D'2 



Tx ~\~ T ntot 

T + T 

J- y r -l ntot 



L {Tx + Tnt °^ 1 (T x + TZ)(Ty+T ntot ) (6) 

where T ntot — T„ F +T nB /Gp and r represents the integra- 
tion time needed to obtain a sky- load measurement. Note 
that the above expressions are relative to power output 
and sensitivity of the complete radiometer, i.e. after aver- 
aging the two radiometer legs; the main advantage of this 
averaging is that any first-order dependency of the knee 
frequency on the mismatch in signal amplitudes after the 
phase switch is cancelled. 

From equation (||) we see that the radiometer sensi- 
tivity is not fully independent of the reference load tem- 
perature, T y , because the signals in the two radiometer 
branches are correlated upstream of the back-end ampli- 
fication. On the other hand for LFI typical parameters, 

is of the order 1(T 4 - 1CT 5 . 



the term 



(T nB /c F y 



Therefore the sensitivity ATo| r = r * in equation (^) can be 
approximated by: 



AT L= r * 



(T,+T ntot ) 



(7) 



which is independent from T y . In the following of this pa- 
per we will always refer to the radiometer sensitivity cal- 
culated for r = Vq even if not explicitly indicated. 

Slightly unbalanced radiometer. We will now consider the 
case of a slight imbalance in the radiometer front-end, 
assuming G Bl = G B , 2 = 1 and T„ Bi = T nB2 = 0. The 
impact of 1/f noise from the back-end amplifiers will be 



analysed in section 3.3 



The first-order approximation of a slightly unbalanced 
front- end can be described mathematically by setting: 



Gf 1 — G 

Tn Fl = 1 

A 1 = l- 

<t> = £0 



G F2 = G(l + c G ) 
i T nF2 = T„(l + e T „) 
A 2 = 1 — e A2 

e 2 = Tr{l + ee 2 ) 



(8) 



In equation (g) we assume that €q < 0.3 (i.e. we as- 
sume a gain mismatch < 1 dB) and that all other e param- 
eters are < 0.1. Then we expand equations (^) and (B.l) 
in series at the second order in eq and at the first order 
in the other parameters, obtaining the following approxi- 
mate equations for the average power output, p, and the 
sensitivity AT (the parameters ct\ through 04 are defined 
in appendix ^) : 



V ~ Po t 1 + V 2 ( £ G - cti£ 2 G - (sa 1 + £A 2 )/2 + a 2 £T„ 

AT « AT (l + 1/4 a 3 e T „ + 1/16 a 4 4) ~ 
w AT (l + l/4a 3 e T J 



(9) 



where we have neglected the term 1/16 04 t 2 G (of the order 
of ~ 10~ 3 e G ) with respect to l/ia^eT,, (of the order of 
- 0.4 e T J. 

From equation (|^) it is apparent that the main ad- 
ditional contribution to the ideal sensitivity comes from 
the unbalance in the front-end noise temperature. For 
LFI typical values the sensitivity is degraded by a fac- 
tor in the range (0.4 — 0.5) x 6t„ which means, in other 
words, that with a noise temperature match better than 
5% it is possible to maintain the sensitivity degradation 
at levels below 1%. A second consideration is that small 
non-idealities in the radiometer chains introduce only a 
very weak dependence of the sensitivity on the reference 
load temperature; in fact from equation (m it follows that 



BAT 



9T U 



T n (T m +T 

2/3r 2(T y +T„) a 



I PS 1.5 X 10" 



3. Susceptibility to various systematic effects 

In this section we study the sensitivity to some potential 
systematic effects arising from fluctuations in the radiome- 
ters; in fact fluctuations in any of the terms appearing in 
equation (|J) will lead to a change in the observed signal 
which can mimic a "true" sky fluctuation. If we denote 
with AT Gq the spurious signal fluctuation induced by a 
variation in a generic radiometer parameter w we have 
that: 



dp A „ dp 



(10) 



In the following sections of this paper, we calculate the 
magnitude of these effects under the assumption that the 
fluctuations in the various parameters are uncorrelated. 
Before calculating these effects we briefly examine the ex- 
pected magnitude of gain and noise temperature fluctua- 
tions. 

3.1. Intrinsic HEMT amplifier noise characteristics 

Cryogenic HEMT amplifiers are well known to have 1/f 
type gain fluctuations, and from this we can infer that 
they also have 1/f type fluctuations in noise temperature 
(Pospieszalski |1989| , Wollack |l995| Jarosik |1996| ). The level 
of these fluctuations can vary considerably among ampli- 
fiers and depends on the details of device fabrication, de- 
vice size, circuit design, and other factors. Because of this, 
we will adopt an empirical model for the fluctuations. We 
can write the 1/f spectrum of the gain fluctuations as: 



(11) 



AG(f) = C 

g V7 

where C represents a constant normalization factor. 
Similarly, we can write the noise temperature fluctuations 
as 

AT n (/) _ A 

Tn VT 

where A is the normalization constant for noise tempera- 
ture fluctuations. Here AT(/)/T has units of Hz -1 / 2 and 



(12) 
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A is dimensionless. From elementary statistical consider- 
ations, we can infer that A — C/2\fW s , where N s is the 
number of stages in the amplifier. 

A normalization of A ~ 1.8 x 10 -5 (see above refer- 
ences) is appropriate for the 30 and 44 GHz radiometers. 
For the 70 and 100 GHz radiometers it will be necessary 
to use HEMT devices with a smaller gate width to achieve 
the lowest amplifier noise figure. We expect that the gate 
widths will be roughly 1/2 that of the devices used for 
the lower frequency radiometers and this will lead to fluc- 
tuations that are roughly a factor of v2 higher. For the 
higher frequencies we will therefore adopt a normaliza- 
tion of A — 2.5 x 10~ 5 . Note that values for A given here 
should be regarded as estimates rather than precise values, 
because in general A will be different for any particular 
device and will generally depend on the physical temper- 
ature of the amplifier. 

While the calculations in this paper assume the above 
simple functional form for AT(/), it is straightforward 
to repeat the calculations with a more detailed spectral 
shape. 

3.2. Front-end amplifier fluctuations 

3.2.1. Sensitivity to noise temperature fluctuations 

In this section, we will calculate the change in the output 
signal for a small change in the noise temperature of the 
front-end amplifiers. Using equation (|l0|) we have that the 
change in input signal induced by a fluctuation in T nF is 
given by: 



AT cq (f) = AT n (f) 



dp/dT n 



dp/dT x 



(13) 



Because both amplifiers (which have uncorrelated 
noise) can contribute to the change in input signal, we 
have that the change in input signal induced by noise tem- 
perature fluctuations in both amplifiers is: 



where AT (/) = y^(T x + T n ). Now with simple algebra 

we can rewrite equation (^) in terms of which equals 
A/yfJ (see equation (jjj)). A further expansion in e yields: 

fk(T„) = / fe °(T„) (1 - 1/2 a 3 e Tn + 1/8 (r + 3) e%) (17) 



A(l-r)T n 
2 (T x +T„) 



and «3 is defined in ap- 



where /°(T») = 
pendix ^|. 

Let us now consider first the ideal case, in which 
fk(T n ) = fk(T n )- Assuming a 20% bandwidth and a ther- 
modynamic sky temperature of 2.7 K, we can calculate 
knee frequencies for several choices of Ty and T n . Results 
are summarised graphically in Figure pi where we show 
contour plots of constant fk as a function of T y (given in 
thermodynamic temperature) and T n . The contour rela- 
tive to fk = /spin = 0.0166 Hz is shown with a dashed 
line. The graphs also show (for each frequency) the range 
of typical noise temperature values (grey area) and the 
nominal reference load temperature (4 K — double hori- 
zontal line) . These results show that at all frequencies the 
expected knee frequency arising from noise temperature 
fluctuations is /°(T„) ^3—5 mHz, i.e. about 4 to 7 
times less than the spin frequency. 

In the limits of our approximations (i.e. a noise tem- 
perature match better than 10% and a gain match better 
than 1 dB) radiometer non idealities determine a correc- 
tion to the zero-order knee frequency that is within ±10%. 

3.2.2. Sensitivity to gain fluctuations 

In this section we calculate the radiometer sensitivity to 
front-end gain fluctuations and show that with proper gain 
modulation the effect is negligible compared to the effect 
induced by fluctuations in T n . 

Proceeding similarly to the previous section the post- 
detection knee frequency of 1// noise caused by gain fluc- 
tuations in the front-end amplifiers (in the case of a per- 
fectly balanced radiometer) is: 



AT eq (f) = AT„(/) 



AT„(/) 



/ dp/dT nFi 
V dp/dT x 



fdp/OTr. 



\ dp/dT x 



V2 



3 + r 
16 



(14) 



If the reference load is at the same temperature as the 
sky, then r = 1 and the effect vanishes. 

Now we calculate the post-detection frequency, fk, at 
which the contributions from noise temperature fluctua- 
tions are equal to the white noise from an ideal radiometer: 



AT eq (/ fc )| r=r » = AT(f) 



(15) 



From equation (j^) we can write the sensitivity spectral 
density as: 



AT(/)=ATo(/)(1 + 1/4o3 6tJ 



(16) 



f° k (G) 



T n - r(T y + T n )Y 



T n f 



(T v + T n y 



(18) 



which shows that with proper gain modulation (r = r$) 
the radiometer is insensitive to gain fluctuations, i.e. 
fk(G) = 0. If we calculate fk(G) in the general case with 
the usual series expansion we have: 



/ fc (G)l , « 



C(T x -T v )(2T n + T x +T y ) 



n 2 



; (T n + T x ) {T n + T y ) 



(19) 



From equation (19) we see that for typical LFI pa- 
rameters and cq ~ ±0.15 (which corresponds to a gain 
mismatch of the order of ~ ± 0.5 dB) the knee frequency 
fk(G)\r=r* < 1 mHz, which indicates that with proper 
gain modulation the 1// noise is still dominated by noise 
temperature fluctuations even if the radiometer is slightly 
unbalanced. 
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Fig. 3. curves of equal fk (in Hz) on the plane T y (K, thermodynamic temperature), T n assuming a thermodynamic 
sky temperature of 2.7 K. Each panel refers to a different frequency channel. The dashed contour refers to values for 
which the knee frequency is equal to the spin frequency (f sp i n = 0.0166 Hz). The graphs also show the range of typical 
LFI noise temperature values (grey area) and the nominal reference load temperature (4 K — double horizontal line) . 



In LFI baseline we foresee a software implementation 
of the gain modulation factor r, so that the sky and ref- 
erence load signals will be detected, converted to digital 
and then subtracted after multiplication of the reference 
load temperature by r. 

Although r will not suffer any fluctuations, it will be in 
general different, at any time, from the ideal value r* . This 
has an impact on the radiometer 1// noise at two different 
levels: (i) /£(G) will be different from zero (because the 
effect of gain oscillations will not be cancelled completely) 
and (ii) f®(T n ) will decrease or increase depending on the 
sign of e r (because r will be closer or farther from 1 ) . 

In order to evaluate the effect of a slight deviation of 
r from r* let us consider a reference value of r that nulls 
or makes very close to zero the output at a certain instant 
to (i.e. r(to) = r*(to)); therefore at a generic time t the 



parameter r(t) can be written as: r(t) = r*(io)(l + £r{t)), 
where e r [t) <C 1. 

Considering that for r ~ r* we have /"(G) ~ 

1 g 2^ [T.+T^-r(T y +T„)] 2 ^ equation we obtain the 

following expressions for fk{T n ) and fk{G): 



fk(Tn)=0 



AT n 



1 2 



2 (T„ + T x ) 



1 - 



(1 + e r ) {Tn + T x ) 
T n + T„ 



(20) 



If we solve for e r the equation f®(T n ) — f®(G) we find 
two solutions, e ri and e r2 (with e ri < e r2 ) given by: 



Tn (Ty T x 
(T n + T x ) 



[(C/A-l) T n + (C/A)T v y 
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(21) 



[{C/A + 1) T n + (C/A) T y 



which define an interval [e^e^-J such that fk{G) <C 
f k (T n ) for e ri < e r < e r2 and f k (G) > f k {T n ) otherwise. 



From equation (21) it is apparent that the width of the in- 



terval [e ri ,e r2 ] is smaller for the high frequency channels, 
characterised by higher values of the noise temperature; 
therefore the requirement on the gain modulation factor 
accuracy is determined by the 100 GHz channel. 

Figure ^ shows the behaviour of e ri and e r2 versus 
the ratio between gain and noise temperature fluctuations 
(i.e. C/A = 2 y/Ns) for the 100 GHz channel. The region 
between the two curves corresponds to values of e r for 
which the noise temperature 1// fluctuations dominate 
over gain instabilities. From the figure it is apparent that 
1// noise arising from gain fluctuations do not dominate 
if the gain modulation factor is accurate at the level of 
< ±0.2%. 

Let us estimate the uncertainty in r introduced by the 
largest expected fluctuation in the signal, i.e. the CMB 
dipole. If the fluctuation induced by the Dipole were larger 
than the requirement, then it would be necessary to cal- 
culate r almost in real-time, in sub-minute time scales. 
From a simple calculation of Sr/r = ST<np \ e / (T x +T n ) 
it follows that 5r/r ~ 10~ 4 , which indicates that it will 
be possible to update r on longer time-scales (of the or- 
der of few days) to account for variations caused by slow 
instrumental drifts. 



T, 



+ T n (l + ef-r T y + T n (l-e) 2 



} (23) 



The change in output caused by a small change in gain 
with the presence of such non-ideal isolation is: 

Op 
dG 

Now, with the correct choice of r, we can make the 
sensitivity to gain fluctuations vanish and null the average 
output of the radiometer simultaneously: 

r = T - + r " (1 + £) : (24) 



T y + T n (l-eY 

Note that this freedom of choice for r is not available to 
correlation radiometer designs without gain modulation, 
and therefore they are subject to 1// fluctuations on the 
non-isolated fraction of the signal. 

3.3. Sensitivity to back-end amplifier gain fluctuations 

In this section we analyse the radiometer sensitivity to 
gain fluctuations of the back end amplifiers. Considering 
that the back-end noise is suppressed by a factor ~ 10 3 
by the front-end amplification we can neglect the contri- 
bution of the noise temperature fluctuations. 

Let us consider the case of a perfectly balanced front- 
end and suppose that the gain value in the two states of 
the phase switch is different by AGs- 

The power output from a single radiometer leg will be 
then given by: 



P = Px -rp y 



0.02 



0.01 



-0.01 



-0.02 




C/A 



Fig. 4. Behaviour of e ri and e r2 (calculated from equa- 
tion (21)) as a function of C/A for the 100 GHz channel. 
From the figure if follows that fk(G) < fk(T n ) (i.e. noise 
temperature fluctuations dominate) if the gain modula- 
tion factor is stable at the level of < ± 0.2 %. 

Let us now analyse the effect of a small leakage across 
the first hybrid coupler, in the balanced radiometer ap- 
proximation; if e is the fraction of a signal applied to one 
port that shows up in the isolated port, then we can write 
the averaged power output, p, as: 



p = akl3G\T x + T n {l 



T, 



+ T„(l-e) 2 ]} (22) 



p x = ak(3G F G B (T x + T np +T nB /G F ) (25) 



p y = akf3G F (G B + AG B ) (T a + T np + T n JG F ) 

A straightforward analysis shows that the fluctuation 
AGb will cause an equivalent spurious signal 

AT eq = -r {T y + T nF + T n JG F ) 



Gi 



(26) 



Substituting AGb/Gb — C/^/J, r — Tq and setting 
Ar c q U i v equal to the single radiometer leg spectral density 
(equal to 2f3~ 1 / 2 (T x + T ntot ), see equation (0)) we can 
calculate the post-detection knee-frequency as: 



(27) 



f k {G B ) = pN s A 2 



Using typical values for /?, N s and A in equation ( |27[ ) 
it follows that the post-detection knee frequency arising 
from back-end gain fluctuations can be much greater than 
1 Hz. Therefore to eliminate the contribution of the back- 
end 1/ / noise we need a phase switch frequency such that 
/sw 3> /fe(Gs), so that the gain will be constant during 
the post-detection integration time. 

For the LFI radiometers we have chosen a phase switch 
frequency of ~ 4 KHz which leads to a small back-end 1// 
contribution to the overall radiometer knee frequency for 
expected values of back-end gain fluctuations. 
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3.4. Sensitivity to reference load fluctuations 

Proceeding as in the previous section we have that the 
signal mimicked by reference load fluctuations, AT y , is 
given by AT eq = -rAT y . 

If we consider random fluctuations in the reference load 
signal with a spectral density given by AT y (f) then we 
have that these fluctuations equal the white noise for: 



AT y (f) = AT°(/) 



1 



"3 CT„ 



r) 16 



(28) 



T n ) and 0:3 is defined in 



where AT°(f) = 
equation ( |C.l| ). 

This provides an upper limit to allowed random fluc- 
tuations in T y in order not to dominate the noise of 
the radiometer. For typical LFI parameters, this upper 
limit is of the order of AT y ~ 1.8 x lO^K/VHz" at 
30 GHz, slightly increasing for the higher frequency chan- 
nels (3.2 x lO^K/VHz" at 100 GHz). If we now consider 
the reference load in a global systematic error budget for 
LFI (the details of which are outside the scope of our pa- 
per) this implies requirements that are about one order of 
magnitude smaller than the above values. 

If we consider non random fluctuations in the reference 
signal then the error induced by these variations must be 
much less than the sensitivity per pixel in the final maps. 
This error, that depends on the spectral behaviour of the 
spurious signal, on the satellite scanning strategy and on 
the data analysis procedures to build the sky maps, can be 
estima ted using the approach described in Mcnnclla ct al. 
( 2002 ) . If we follow this approach considering that the pa- 
rameter r is of order unity, we have that spin synchronous 
variations of the reference signal will be transferred di- 
rectly to the final maps, while slower fluctuations will be 
damped by a factor of the order of 10~ 2 - 10~ 3 by the 
measurement strategy and by data analysis. 

4. Conclusions 

In this paper we have discussed the pseudo-correlation ar- 
chitecture adopted for the radiometers of the Planck-LFI 
instrument and we have studied the sensitivity of the mea- 
sured signal to various systematic effects. In our treatment 
we have considered both the ideal case of a perfectly bal- 
anced radiometer and the effect of small mismatches in 
the various radiometer parameters. 

The first result is that the radiometer sensitivity does 
not depend on the level of the reference load temperature; 
even in the case of a slight imbalance in the radiome- 
ter parameters the dependence on T y is at the level of 
dAT/dTy ~ 10 -5 which is negligible. The only mismatch 
which has a first-order impact on AT is the noise temper- 
ature mismatch of the front-end amplifiers; our analysis 
shows that it is possible to maintain the sensitivity degra- 
dation below 1% with a noise temperature match better 
than 5%. 

With proper gain modulation (r = rg) the 1/f noise 
in the radiometer output is determined mainly by noise 



temperature fluctuations in the front-end amplifiers, with 
a knee frequency of few mHz, provided that the front-end 
amplifier amplitude match is better that ~ ±0.5 dB. Such 
a high level of 1/f noise suppression depends on the gain 
modulation factor, which must be determined with an ac- 
curacy better than ±0.2%. If the accuracy on r is less than 
±1% then gain fluctuations become the major source of 
1/f noise; if e r is kept in the range ±1% we expect val- 
ues of the knee frequency of the order of 50 mHz which 
can be easily handled by destriping algorithms. The pres- 
ence of a small amount of leakage in the first hybrid does 
not significantly modify this conclusion; with the correct 
choice of r, one can make the sensitivity to gain fluctua- 
tions vanish and null the average output of the radiometer 
simultaneously. 

The effect of gain fluctuations in the back-end ampli- 
fiers can be made negligible by the fast front-end switching 
between sky and reference load signals. The LFI baseline 
of 4096 Hz for the phase switch frequency has been cho- 
sen to guarantee a high suppression level of the 1/f noise 
from back-end amplifiers. 

In general our analysis demonstrates the effectiveness 
of the gain modulation concept applied to this form of 
radiometer. The estimate of the knee frequency given by 
equation (|l7]) is quite low and relatively immune to small 
imperfections in radiometer balance. The modified corre- 
lation radiometer scheme reduces the knee frequency by 
more than two orders of magnitude, compared to a total 
power radiometer of similar bandwidth and intrinsic tran- 
sistor fluctuations. For such small residual knee frequency 
(of ~ 0.1 Hz) it will be possible to remove efficiently the 
effects during data analysis. 

We have also studied the sensitivity of the radimc- 
ter to changes in the reference load temperature, T y . Our 
analysis provides a framework in which thermal stability 
requirements on the LFI reference loads can be evaluated. 

A refinement of the present analysis for the determi- 
nation of fk will be pursued in the future by software sim- 
ulations of the radiometer functions to accurately study 
the combined effect of all components. Finally, laboratory 
measurements of a prototype radiometer working under 
conditions close to those of Planck mission constitute the 
most important checks for the ultimately understanding of 
the behaviour of Planck LFI radiometers regarding the 
1// type noise and possible further effects. Results of pre- 
liminary laboratory measurements performed of Planck- 
LFI prototype radiometers will be presented in forthcom- 
ing publications. 
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Appendix A: Power output 

The radiometer power output at each of the radiometer output legs is defined by: 



Pi (*) = T 9%i { ( 2 n -Bi cos i ( t ) B l ) + g Fl (V2n Fl + x + y) cos (<j> Fl + 4> Bl 
+ \fA[g F2 {\/2 np 2 + x - y) cos (Ox + <\>p 2 + 4>bS) + 
+ (2 n Bl sin (4> Bl ) + g Fl (V^n Fl + x + y) sin (<j> Fl + 4> Bl ) + 

+ \/~AlgF 2 {V2n F2 + x - y) sin (6>i + <j) F2 + 4> Bl )j 4 
-r (2n Bl cos {4> Bl ) + g Fl (V2n Fl + x + y) cos (<f> Fl + 4> Bl )+ 

+ \fM gF 2 (V2nF 2 + x - y) cos (9 2 + 4>f 2 + ^bJ^ + 
+ (2 n Bl sin (4> Bl ) + g Fl {V2n Fl + x + y) sin (cf> Fl + 4> Bl ) + 

+ \^A^g F2 (V2n F2 + x- y)sm(9 2 + <pF 2 + fe)) 



(A.l) 



where a is the constant of proportionality of the square-law detectors. A similar relationship holds for P2(t). 



Appendix B: Sensitivity 

B.l. Outline for sensitivity calculation 



In this appendix we outline the procedure to calculate the LFI radiometer sensitivity (see equation (B.l)). We start from 
the expression of the power output p(t), which is given by equation (A.l) and then we calculate the autocorrelation 
of p(t), denoted as iP p (t) = p(t)p(t + r). Note that we assume that x, y, m, and n% are all uncorrelated gaussian 
variables so that xy = xy and that they all have zero mean. The next step is the calculation of the Fourier transform 
of tpp( T ) (indicated by Q p (f), which can be split in a constant part, Q p (0), and in a fluctuating part Q' p (f). Then we 

calculate the rms voltage w lms = w(t) 2 ^ where w(t) 2 = df H (f)Q' p (f) and H(f) is a rectangular lowpass filter 
of height Hq and width b. The last step is to calculate the change in power output Ap determined by a change in the 



sky signal equal to AT; the sensitivity AT is then calculated by solving the equation w rl 



#o 1/2 Ap. 



B.2. Sensitivity analytical form 

The radiometer sensitivity AT (calculated according to the procedure outlined above) has the following form: 

AT = - [AT? 
2 L 1 



AT 2 2 1 3 



(B.l) 



where ATi and AT2 (which represent the minimum detectable rms signal on each single radiometer leg) are defined 
by the following equations: 



ATi = 
AT 2 = 



~T(A-Br + Cr 2 Y 
2/^ G-rl 
~T(D -Er + Fr 2 y 



2(3t G-rl 
The coefficients A through F have the following form: 
A = a{F'l + a 2 Ty + a 3 T x + a 4 T y + a 5 T x T y + a 6 

with similar relationships for B, C, D, E and F. Coefficients a,j, bj, etc. are defined as follows: 



(B.2) 



(B.3) 
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B.2.1. Definitions of a,j coefficients 



ai 
a 3 



2G 2 
2 1 2 

0-31 T nFl + a 32 T nF2 + a 33 T nBi 
a 3i = 2 



\/G Fl + V A 1 G F 2 COS(01 - (j) Fl + (j) F2 



1 



a.32 — 2 Ax ^2 



^Gf 1 cos(^i - Fl + (^) F2 ) 



033 



cos(0 Fl ) + ^/ Ai Gf 2 cos(6»i + <j> F2 ) 



a.4 



0.41 T nFl + d42 T n2 + d 43 T n 

an = 2 G F 

1 . _ 

a 42 — 2 1 F2 



a/G Fi - y/Ax G F2 cos(0i - Fl + 4>p 2 ) 
\[A~x Gf 2 - \/G~f\~ cos(6»i - </> Fl + 4>f 2 ) 



043 



7 G Fl cos(0 Fl ) - a/ G F2 cos(6»i + </> F2 ) 



a 5 
a 6 



- (G Fl - Ai G F2 ) 2 

- G Fl T„ Fi 



G Fl T , „ Fi + 4cos(</> Fl ) T„ 



+ 2 ^i G F2 T„ F2 



2 cos(0i - <£ Fl + 4>f 2 ) G Fl T nFi + Ai G F2 T n2 + 4 cos(9 1 + 4>f 2 ) T nBi 



+ 2T 2 

1 ^ ^ 71 



B.2.2. Definitions of bj coefficients 



h= AI G 
b 2 = AIG 

b 3 = hi T nFl + b 32 T nF2 + b 33 T nBi 



b 3 i = G Fl y / G Fl + \J A x Gp 2 cos(6»i - <p Fl + 4>p 2 ) \/G Fl + V 'A 2 Gp 2 cos(0 2 - 4>f 1 + 4>f 2 

b 32 = \J Ax A 2 Gf 2 \J Ai Gf 2 + \/G Fl cos(0i - <p Fl + <Pf 2 ) \J A 2 G F2 + \/G Fl cos(0 2 - cp Fl + 4>f 2 

b 33 = 2 v / G Fl cos(</> Fl ) + y/ Ax Gp 2 cos(0i + 4>f 2 ) \/G Fl cos(0 Fl ) + x/ A 2 Gf 2 cos(6> 2 + 4>f 2 ) 
64 = &4i T nFi + 642 T nF2 + 643 T nBi 

&41 = G Fl \/G Fl - \J A x G F2 cos(6»i - 4> Fl + 4>f 2 ) \/G Fl - \J A 2 Gp 2 cos(9 2 - 4> Fl + 4>f, 



b i2 = x/Ax A 2 Gp 2 xj Ax Gp 2 - \]~G F ~ l cos(dx - </> Fl + <Pf 2 ) \J A 2 Gp 2 - \]~G F ~ 1 cos(0 2 - <f> Fl + 4>f 2 ) 



b 43 = 2 
1 



G Fl cos(0 Fl ) - \J Ax G F2 008(01 + <Pf 2 ) \/G Fl cos(0 Fl ) - \J A 2 Gf 2 cos(0 2 + 4>f 2 ) 
b 5 = - (G Fl - Ax G F2 ) (G Fl - A 2 G F2 ) 
b 6 = G Fl T nFi G Fl T nFi + 4 cos(> Fl ) 



+ 



+ G F2 \/Ax A 2 T nF2 yjAx A 2 G F2 T„ F2 + 4T„ Bi cos(0i + F J cos(0 2 + <Pf 2 ) + 
+ 2 G Fl T„ Fi cos(0i - Fl + </> F2 ) cos(0 2 - 4> Fl + <f>F 2 )] + 4 T 2 Bi 



B.2.3. Definitions of c, coefficients 



Cl 
C2 



2/ 2 
2G 2 



12 

C3 



c 4 
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C31 T„ Fi + c 32 T nF2 + c 33 T„ Bi 

1 r 12 

- G Fl a/ G Fl + V ^2 Gp 2 COs(6» 2 - 0Fi + (f)F 2 ) 



C31 



C32 — 2 ^2 Gf 2 



C33 



a/A 2 Gf 2 + \/Gf\ COs(0 2 - 0Fi + 0F 2 

2 



/ Gf7 COS(0 F J + y^4 2 Gf 2 COs(0 2 + 4>F 2 ) 



C41 ?n Fi + C 42 T„p 2 + C43 T„ Bi 
1 



C41 



G 



c 42 — ^ ^ 2 



G Fl - V A 2 G F2 C0S(6» 2 - </>Fi + </>F 2 ) 



\/ A 2 Gf 2 - x/G Fl COs(# 2 - <M + ^F 2 



C43 



/ G Fl COs((fi Fl ) - y^4 2 Gf 2 COs(# 2 + </>F 2 ) 



C5 
C6 



1 2 

- (G Fl - A 2 Gf 2 ) 
1 



— 7: G F , t„ 



2 J Fl J «Fi 

+ 2 



Gfj T„ Fi +4cos(0fJ 



+ 



i A 2 Gf 2 T nF2 [2 cos(fl 2 - Fl + <f> F2 ) 2 G Fl T„ Fi + A 2 Gf 2 T„ F2 + 4 cos(# 2 + </>f 2 ) 2 T„ Bi 



B.2.4. Definitions of coefficients 



di 
d 2 
d 3 



d 5 = 
d 6 = 



2G 2 
2I 2 

d3,l Tn Fl + dz.2 T np2 + ^3,3 T n 

1 



fife,! = 2 G Fi V G Fi - V A 2 G F 2 C0S(6» 2 - 4>F 1 + §F 2 ) 



^3,2 — ^ ^2 Gf 2 



V^GfT - ^Gfj C0S(# 2 - 0Fi + 0F 2 



*3,3 



/ Gfi COs(0fJ - \A4 2 Gf 2 COs(0 2 + </>F 2 ) 



c^4 — £^4,1 T nFi + C?4 :2 T„ F2 + ^4,3 T„ B2 



^4,1 = 2 + V^GfT C0S(6» 2 - 0Fj + </>F 2 ) 



^4,2 — 2 ^ 2 



V^GfT + COs((9 2 - (/> Fl + 0F 2 

1 2 



/ Gf7 COs(0fJ + V A 2 Gp 2 COs(6» 2 + <f>F 2 ) 



di,3 = 

- {G Fl - A 2 Gf 2 ) 2 

1 r 2 

- G Fl T„ Fi Gfi T nFi + 4 cos(<£ Fl ) 



+ 



+ - A 2 Gf 2 T„ F2 



2 cos(fl 2 - <fo + 0f 2 ) 2 Gfi T„ Fi + A 2 G F2 T nF2 + 4cos(0 2 + <fe) 2 T„ B2 



+ 2T„ 2 



B.2.5. Definitions of e, coefficients 



ei = 4/G 
e 2 = 4/G 

63 = ^3,1 T nFi + d 3i2 T„ F2 + d 3:3 T nB2 



e3,i = G Fl vGfi - V^iGfT cos(6»i - <j> Fl + <j)p 2 ) y/G Fl - \J A 2 Gf 2 cos(fl 2 - 4>f-_ 



e3,2 = \JMA 2 G F2 y/ Ax Gf 2 - ^GfT cos(^i - 0f x + </>f 2 ) V ^2 Gf 2 - V^f\ cos(0 2 
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e 4 



63,3 = 2 \/G Fl cos(0 Fl ) - y/A~i G F2 cos(6»i + <j) F2 ) \jG Fl cos(<f> Fl ) - \/A 2 G F . 2 cos(0 2 + §p 2 ) 
e4,i T nFi + e^2 T UF2 + T„ B2 



e4,i = G Fl \jG Fl + y/Ax G F2 cos(#i - <j> Fl + <pp 2 ) \jG Fl + V ' M G F2 coa(0 2 - <p Fl + <Pf 2 
e4,2 = \[A\ A 2 Gp 2 sj A x Gp 2 + \jG Fl cos(9 1 - 4> Fl + 4>p 2 ) y/ A 2 G F . 2 + y/G Fl cos(6 2 - <f> Fl + 4>p 2 ) 



e4,3 = 2 \jG Fl cos(<^Fi) + \/ Ai Gf 2 cos(#i + (f>p 2 ) \jG Fl cos(<p Fl ) + \J A 2 Gp 2 cos(8 2 + 4>p 2 ) 



e 5 — — (G Fl — A\ Gf 2 ) (G Fl — A 2 Gp 2 



+ 



e 6 = G Fl T nFi G Fl T nFi + 4cos((/) Fl ) 

+ Gp 2 y/Ai A 2 T nF2 \y/Ai A 2 G F2 T nF2 + 4T„ B2 008(0! + cj> F2 ) cos(6> 2 + <f> Fa )+ 
+ 2 G Fl T nFi cos(6»i - <f> Fl + <Pp 2 ) cos(e» 2 - 4> Fl + <\>f 2 )] + 4 



B.2.6. Definitions of coefficients 



h= 2 1 2 

h= 2G 2 

fa, = h.i T nFi + /3,2 T nF2 + f 3j3 T, 

1 



/3,i - 2 Gf x 
h,2 — 2 M Gp 2 



G Fl - y/ Ai G F2 cos(^i - <f> Fl + <j)p 2 ) 



\J Ai Gf 2 ~ \JGf-i cos(6»i - <f> Fl + (j>p 2 ) 
/3,3 = cos(4> Fl ) - xj Ai Gp 2 cos(0i + 4>F 2 ) 

fi = fi.l T nFl + fi,2 Tn F2 + /4,3 T nB2 

h,i = 7; G Fl Vg^ + \/Ai G F2 cos(6>i - 4> Fl + 4>f 2 ) 

1 r i2 

fi,2 = 2 A i g f 2 \yAi G F2 + \JG Fl cos(#i - 4> Fl + 4>f 2 ) 

l 2 



.A, 3 = COS^fJ + COS(01 + </>f 2 

1 _ 2 

/s = ^ (Gfi - Ai Gf 2 ) 



k = 



1 



Gf I? 



G Fl T„ Fi +4cos(^ 1 ) 2 T„ 



+ 2 ^i Gf 2 r„ F2 



2cos(^i - <f> Fl + 4>F 2 f G Fl T nFi + Ai G F2 T nF2 +4008(0! + <\> F2 f T n 



+ 2T„ 2 
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Appendix C: Definition of a parameters 

The parameters a.\ through a<± used in equations 
( |l7| ) and are defined as: 



(l + r)(T x -T y ) 

OCl 



a 2 



<*3 = T n 



8[T x + T n -r (T n + T n )] 

(1 - r) T n 
T x +T n -r [T n + T n ) 
1 1 



T% T T n Ty -\- T n 



Tx Ty , n ^ 

«4 = — =r (C.l) 

T 4- T 
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